Transcriptional drift in aging cells: A global de-controller.

As cells age, they undergo a remarkable global change: In transcriptional drift, hundreds of genes become overexpressed while hundreds of others become underexpressed. Using archetype modeling and Gene Ontology analysis on data from aging Caenorhabditis elegans worms, we find that the upregulated genes code for sensory proteins upstream of stress responses and downregulated genes are growth- and metabolism-related. We propose a simple mechanistic model for how such global coordination of multi-protein expression levels may be achieved by the binding of a single ligand that concentrates with age. A key implication is that a cell's own responses are part of its aging process, so unlike for wear-and-tear processes, intervention might be able to modulate these effects.

Adaptive DNA amplification of synthetic gene circuit opens a way to overcome cancer chemoresistance.

Drug resistance continues to impede the success of cancer treatments, creating a need for experimental model systems that are broad, yet simple, to allow the identification of mechanisms and novel countermeasures applicable to many cancer types. To address these needs, we investigated a set of engineered mammalian cell lines with synthetic gene circuits integrated into their genome that evolved resistance to Puromycin. We identified DNA amplification as the mechanism underlying drug resistance in 4 out of 6 replicate populations. Triplex-forming oligonucleotide (TFO) treatment combined with Puromycin could efficiently suppress the growth of cell populations with DNA amplification. Similar observations in human cancer cell lines suggest that TFOs could be broadly applicable to mitigate drug resistance, one of the major difficulties in treating cancer.

Cusp bifurcation in a metastatic regulatory network.

Understanding the potential for cancers to metastasize is still relatively unknown. While many predictive methods may use deep learning or stochastic processes, we highlight a long standing mathematical concept that may be useful for modeling metastatic breast cancer systems. Ordinary differential equations (ODEs) can model cell state transitions by considering the pertinent environmental variables as well as the paths systems take over time. Bifurcation theory is a branch of dynamical systems which studies changes in the behavior of an ODE system while one or more parameters are varied. Many studies have applied concepts in one-parameter bifurcation theory to model biological network dynamics, and cell division. However, studies of two-parameter bifurcations are much more rare. Two-parameter bifurcations have not been studied in metastatic systems. Here we show how a specific two-parameter bifurcation phenomenon called a cusp bifurcation separates two qualitatively different metastatic cell state transitions modalities and propose a new perspective on defining such transitions based on mathematical theory. We hope the observations and verification methods detailed here may help in the understanding of metastatic potential from a basic biological perspective and in clinical settings.

Nonmonotone invasion landscape by noise-aware control of metastasis activator levels.

A major pharmacological assumption is that lowering disease-promoting protein levels is generally beneficial. For example, inhibiting metastasis activator BACH1 is proposed to decrease cancer metastases. Testing such assumptions requires approaches to measure disease phenotypes while precisely adjusting disease-promoting protein levels. Here we developed a two-step strategy to integrate protein-level tuning, noise-aware synthetic gene circuits into a well-defined human genomic safe harbor locus. Unexpectedly, engineered MDA-MB-231 metastatic human breast cancer cells become more, then less and then more invasive as we tune BACH1 levels up, irrespective of the native BACH1. BACH1 expression shifts in invading cells, and expression of BACH1's transcriptional targets confirm BACH1's nonmonotone phenotypic and regulatory effects. Thus, chemical inhibition of BACH1 could have unwanted effects on invasion. Additionally, BACH1's expression variability aids invasion at high BACH1 expression. Overall, precisely engineered, noise-aware protein-level control is necessary and important to unravel disease effects of genes to improve clinical drug efficacy.

MaxCal can infer models from coupled stochastic trajectories of gene expression and cell division.

Gene expression is inherently noisy due to small numbers of proteins and nucleic acids inside a cell. Likewise, cell division is stochastic, particularly when tracking at the level of a single cell. The two can be coupled when gene expression affects the rate of cell division. Single-cell time-lapse experiments can measure both fluctuations by simultaneously recording protein levels inside a cell and its stochastic division. These information-rich noisy trajectory data sets can be harnessed to learn about the underlying molecular and cellular details that are often not known a priori. A critical question is: How can we infer a model given data where fluctuations at two levels-gene expression and cell division-are intricately convoluted? We show the principle of maximum caliber (MaxCal)-integrated within a Bayesian framework-can be used to infer several cellular and molecular details (division rates, protein production, and degradation rates) from these coupled stochastic trajectories (CSTs). We demonstrate this proof of concept using synthetic data generated from a known model. An additional challenge in data analysis is that trajectories are often not in protein numbers, but in noisy fluorescence that depends on protein number in a probabilistic manner. We again show that MaxCal can infer important molecular and cellular rates even when data are in fluorescence, another example of CST with three confounding factors-gene expression noise, cell division noise, and fluorescence distortion-all coupled. Our approach will provide guidance to build models in synthetic biology experiments as well as general biological systems where examples of CSTs are abundant.

Immune mechanisms shape the clonal landscape during early progression of prostate cancer.

Understanding the role of the immune microenvironment in modulating intratumor heterogeneity is essential for effective cancer therapies. Using multicolor lineage tracing in genetically engineered mouse models and single-cell transcriptomics, we show that slowly progressing tumors contain a multiclonal landscape of relatively homogeneous subpopulations within a well-organized tumor microenvironment. In more advanced and aggressive tumors, however, the multiclonal landscape develops into competing dominant and minor clones accompanied by a disordered microenvironment. We demonstrate that this dominant/minor landscape is associated with differential immunoediting, in which minor clones are marked by an increased expression of IFNγ-response genes and the T cell-activating chemokines Cxcl9 and Cxcl11. Furthermore, immunomodulation of the IFNγ pathway can rescue minor clones from elimination. Notably, the immune-specific gene signature of minor clones exhibits a prognostic value for biochemical recurrence-free survival in human prostate cancer. These findings suggest new immunotherapy approaches for modulating clonal fitness and tumor progression in prostate cancer.

A luminal intermediate cell state maintains long-term prostate homeostasis and contributes to tumorigenesis.

Cellular heterogeneity poses tremendous challenges for developing cell-targeted therapies and biomarkers of clinically significant prostate cancer. The origins of this heterogeneity within normal adult and aging tissue remain unknown, leaving cellular states and transcriptional programs that allow expansions of malignant clones unidentified. To define cell states that contribute to early cancer development, we performed clonal analyses and single cell transcriptomics of normal prostate from genetically-engineered mouse models. We uncovered a luminal transcriptional state with a unique "basal-like" Wnt/p63 signaling ( luminal intermediate , LumI) which contributes to the maintenance of long-term prostate homeostasis. Moreover, LumI cells greatly expand during early stages of tumorigenesis in several mouse models of prostate cancer. Genetic ablation of p63 in vivo in luminal cells reduced the formation of aggressive clones in mouse prostate tumor models. Finally, the LumI cells and Wnt signaling appear to significantly increase in human aging prostate and prostate cancer samples, highlighting the importance of this hybrid cell state for human pathologies with potential translational impact.

New antivirals exploit viral feedback tricks for a cure without resistance.

Many approved drugs, including antivirals, are small-molecule inhibitors of disease-causing proteins. Such inhibitors often elicit resistance during treatment. Chaturvedi et al. propose new, feedback-disruptor (FD) antivirals that efficiently cure infected cells from viruses and minimize the chance of resistance, providing a new paradigm to treat viral infections and possibly other diseases.

Drug-dependent growth curve reshaping reveals mechanisms of antifungal resistance in Saccharomyces cerevisiae.

Microbial drug resistance is an emerging global challenge. Current drug resistance assays tend to be simplistic, ignoring complexities of resistance manifestations and mechanisms, such as multicellularity. Here, we characterize multicellular and molecular sources of drug resistance upon deleting the AMN1 gene responsible for clumping multicellularity in a budding yeast strain, causing it to become unicellular. Computational analysis of growth curve changes upon drug treatment indicates that the unicellular strain is more sensitive to four common antifungals. Quantitative models uncover entwined multicellular and molecular processes underlying these differences in sensitivity and suggest AMN1 as an antifungal target in clumping pathogenic yeasts. Similar experimental and mathematical modeling pipelines could reveal multicellular and molecular drug resistance mechanisms, leading to more effective treatments against various microbial infections and possibly even cancers.

Gene-circuit therapy on the horizon: synthetic biology tools for engineered therapeutics.

Therapeutic genome modification requires precise control over the introduced therapeutic functions. Current approaches of gene and cell therapy fail to deliver such command and rely on semi-quantitative methods with limited influence on timing, contextuality and levels of transgene expression, and hence on therapeutic function. Synthetic biology offers new opportunities for quantitative functionality in designing therapeutic systems and their components. Here, we discuss synthetic biology tools in their therapeutic context, with examples of proof-of-principle and clinical applications of engineered synthetic biomolecules and higher-order functional systems, i.e. gene circuits. We also present the prospects of future development towards advanced gene-circuit therapy.

Reliably Engineering and Controlling Stable Optogenetic Gene Circuits in Mammalian Cells.

Reliable gene expression control in mammalian cells requires tools with high fold change, low noise, and determined input-to-output transfer functions, regardless of the method used. Toward this goal, optogenetic gene expression systems have gained much attention over the past decade for spatiotemporal control of protein levels in mammalian cells. However, most existing circuits controlling light-induced gene expression vary in architecture, are expressed from plasmids, and utilize variable optogenetic equipment, creating a need to explore characterization and standardization of optogenetic components in stable cell lines. Here, the study provides an experimental pipeline of reliable gene circuit construction, integration, and characterization for controlling light-inducible gene expression in mammalian cells, using a negative feedback optogenetic circuit as a case example. The protocols also illustrate how standardizing optogenetic equipment and light regimes can reliably reveal gene circuit features such as gene expression noise and protein expression magnitude. Lastly, this paper may be of use for laboratories unfamiliar with optogenetics who wish to adopt such technology. The pipeline described here should apply for other optogenetic circuits in mammalian cells, allowing for more reliable, detailed characterization and control of gene expression at the transcriptional, proteomic, and ultimately phenotypic level in mammalian cells.

From Bench to Keyboard and Back Again: A Brief History of Lambda Phage Modeling.

Cellular decision making is the process whereby cells choose one developmental pathway from multiple possible ones, either spontaneously or due to environmental stimuli. Examples in various cell types suggest an almost inexhaustible plethora of underlying molecular mechanisms. In general, cellular decisions rely on the gene regulatory network, which integrates external signals to drive cell fate choice. The search for general principles of such a process benefits from appropriate biological model systems that reveal how and why certain gene regulatory mechanisms drive specific cellular decisions according to ecological context and evolutionary outcomes. In this article, we review the historical and ongoing development of the phage lambda lysis-lysogeny decision as a model system to investigate all aspects of cellular decision making. The unique generality, simplicity, and richness of phage lambda decision making render it a constant source ofmathematical modeling-aided inspiration across all of biology. We discuss the origins and progress of quantitative phage lambda modeling from the 1950s until today, as well as its possible future directions. We provide examples of how modeling enabled methods and theory development, leading to new biological insights by revealing gaps in the theory and pinpointing areas requiring further experimental investigation. Overall, we highlight the utility of theoretical approaches both as predictive tools, to forecast the outcome of novel experiments, and as explanatory tools, to elucidate the natural processes underlying experimental data.

Binary Expression Enhances Reliability of Messaging in Gene Networks.

The promoter state of a gene and its expression levels are modulated by the amounts of transcription factors interacting with its regulatory regions. Hence, one may interpret a gene network as a communicating system in which the state of the promoter of a gene (the source) is communicated by the amounts of transcription factors that it expresses (the message) to modulate the state of the promoter and expression levels of another gene (the receptor). The reliability of the gene network dynamics can be quantified by Shannon's entropy of the message and the mutual information between the message and the promoter state. Here we consider a stochastic model for a binary gene and use its exact steady state solutions to calculate the entropy and mutual information. We show that a slow switching promoter with long and equally standing ON and OFF states maximizes the mutual information and reduces entropy. That is a binary gene expression regime generating a high variance message governed by a bimodal probability distribution with peaks of the same height. Our results indicate that Shannon's theory can be a powerful framework for understanding how bursty gene expression conciliates with the striking spatio-temporal precision exhibited in pattern formation of developing organisms.

How Do Cells Adapt? Stories Told in Landscapes.

Cells adapt to changing environments. Perturb a cell and it returns to a point of homeostasis. Perturb a population and it evolves toward a fitness peak. We review quantitative models of the forces of adaptation and their visualizations on landscapes. While some adaptations result from single mutations or few-gene effects, others are more cooperative, more delocalized in the genome, and more universal and physical. For example, homeostasis and evolution depend on protein folding and aggregation, energy and protein production, protein diffusion, molecular motor speeds and efficiencies, and protein expression levels. Models provide a way to learn about the fitness of cells and cell populations by making and testing hypotheses.

Multiplexed Gene Expression Tuning with Orthogonal Synthetic Gene Circuits.

Understanding the individual and joint contribution of multiple protein levels toward a phenotype requires precise and tunable multigene expression control. Here we introduce a pair of mammalian synthetic gene circuits that linearly and orthogonally control the expression of two reporter genes in mammalian cells with low variability in response to chemical inducers introduced into the growth medium. These gene expression systems can be used to simultaneously probe the individual and joint effects of two gene product concentrations on a cellular phenotype in basic research or biomedical applications.

Evolutionary regain of lost gene circuit function.

Evolutionary reversibility-the ability to regain a lost function-is an important problem both in evolutionary and synthetic biology, where repairing natural or synthetic systems broken by evolutionary processes may be valuable. Here, we use a synthetic positive-feedback (PF) gene circuit integrated into haploid Saccharomyces cerevisiae cells to test if the population can restore lost PF function. In previous evolution experiments, mutations in a gene eliminated the fitness costs of PF activation. Since PF activation also provides drug resistance, exposing such compromised or broken mutants to both drug and inducer should create selection pressure to regain drug resistance and possibly PF function. Indeed, evolving 7 PF mutant strains in the presence of drug revealed 3 adaptation scenarios through genomic, PF-external mutations that elevate PF basal expression, possibly by affecting transcription, translation, degradation, and other fundamental cellular processes. Nonfunctional mutants gained drug resistance without ever developing high expression, while quasifunctional and dysfunctional PF mutants developed high expression nongenetically, which then diminished, although more slowly for dysfunctional mutants where revertant clones arose. These results highlight how intracellular context, such as the growth rate, can affect regulatory network dynamics and evolutionary dynamics, which has important consequences for understanding the evolution of drug resistance and developing future synthetic biology applications.